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Abstract 

The massive continuum limit of the 1 + 1 dimensional 0(2) nonlinear 
(T-model (XY model) is studied using its equivalence to the Sine-Gordon 
model at its asymptotically free point. It is shown that leading lattice 
artifacts are universal but they vanish only as inverse powers of the 
logarithm of the correlation length. Such leading artifacts are calculated 
for the case of the scattering phase shifts and the correlation function 
of the Noether current using the bootstrap S-matrix and perturbation 
theory respectively. 
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1 Introduction 



In this paper we study the properties of the two-dimensional 0(2) nonhnear cr- 
model, better known as the XY model. This model has been the subject of extensive 
theoretical and numerical analysis, starting with the seminal papers of Kosterlitz 
and Thouless (KT) For a review of KT theory, see 

Analitical work is usually based on a series of mappings that starts at the original 
(lattice) XY model and arrives at the Sine-Gordon (SG) model or its fermionic 
equivalent, a (deformed) version of the chiral Gross-Neveu (CGN) model. This 
latter formulation is most useful if one wants to study questions related to the 
dynamically generated SU{2) symmetry of the model. 

Most papers on the XY model study its properties interesting for Statistical 
Physics, in particular the pecularities of the KT phase transition, which is of infinite 
order. In this paper we look at the XY model as an example of 1 -|- 1 dimensional 
relativistic Quantum Field Theory. More precisely, we study the massive continuum 
limit of the lattice theory, which, in the language of Statistical Physics, means that 
we approach the KT phase transition point from the high temperature phase. 

Treating the XY model as the n = 2 member of the family of 0{n) nonlinear 
cr- models gives additional insights, since a lot is known about the n > 3 models 0. 
More importantly, using the SG language, we show that the approach to the contin- 
uum limit in this model is much slower than in most other lattice models. Lattice 
artifacts vanish in this model, instead of the usual Symanzik type behaviour |^ 
(i.e. integer powers of the lattice spacing), as inverse powers of the logarithm of the 
lattice spacing only. On the other hand, we can show that the leading artifacts are 
universal and calculable. Our main result is Eqs. (p^ and (Q) in Section 4, which 
allow us to calculate leading lattice artifacts in terms of SG data. 

We can make use of the fact that the SG model is exactly solvable and its boot- 
strap S-matrix is exactly known. We calculate the leading artifacts for the scattering 
phase shifts using the bootstrap results. An alternative method is perturbation the- 
ory (PT). Since the SG model is asymptotically free if we use suitable expansion 
parameters 0, the methods of renormalization group (RG) improved PT are thus 
available. 

In Section 2 we review the relation of the XY model to the 0{n) models with 
n > 3 and describe the chain of mappings leading from the XY model to the SG 
model and its equivalent fermionic formulation. 

In Section 3 we recall the analysis of the phase diagram of the model in the 
vicinity of the KT phase transition point. This is described in the SG language. 

In Section 4 we explain how to calculate the lattice artifacts and apply this to 
the case of the scattering phase shifts. 

Finally in Section 5 we calculate the lattice artifacts for the two-point correlation 
function of the Noether current corresponding to the 0(2) symmetry. Here we 
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use the method of RG improved PT. To calculate the value of a non-perturbative 
constant needed here, we also consider the system in the presence of an external 
field coupled to the Noether charge. (This calculation is analogous to the one used 
previously to determine the M/A ratio for the 0{n) models 0.) We give here a 
parameter-free two-loop formula for the lattice artifacts. 

A precision MC study of the massive continuum limit of the 0(2) model will be 
described in a forthcoming paper 0. 



2 From the XY model to the Sine-Gordon model 
and beyond 

In this section we describe in some detail the chain of mappings starting with the 
XY model and ending at the SG model and its fermionic equivalent. 

We can treat the XY model as the n = 2 member of the family of 0{n) nonlinear 
cr-models with Lagrangean 

CO{n) ^ _}_QgaQga . 3^8'' = I, a = 1, 2, . . . Tl. (1) 

The n > 3 models are known to be integrable. Polyakov P] and Liischer p have 
shown the existence of respectively local and nonlocal higher spin conserved charges, 
whose existence implies quantum integrability. Assuming the spectrum of the model 
consisted of an 0{n) vector multiplet of massive particles the exact S- matrix of the 
n > 3 models was found by bootstrap methods 0]: 



where 



Sab;cdi^) = (^li(^)^abScd + 0-2i6)6acSbd + Cr3id)SadSbc, (2) 



in -9 {n- 2)9 - 2ni ' 

M^) = (n-2).. („_,;i,^^ . (3) 

s^'^H9\ 

a^{9) = -2m ^' 



{n - 2)9 - 2ni 
and the 'isospin 2' phase shift s*-^-* is given by 



s(2)(^) = -exp|2zj^"^sina;^ir„(cu)| (4) 
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with 



-nuj _^ g-27r — 



1 + e- 



(5) 



Much less is known about the 0(2) model. A simple observation is that @ and 
also (^) have a smooth n —>■ 2 limit. It is natural to assume that the 0(2) model 
is also integrable, its spectrum consists of a single 0(2) doublet of massive particles 
whose scattering is indeed described by the — > 2 limit of the S-matrix (|^). 

Although taking the formal n —>■ 2 limit of the bootstrap results valid for n > 3 
is not convincing in itself, the conclusion turns out to be correct because as we will 
see it also follows from the Kosterlitz-Thouless theory |I| of the XY model. 

Before turning to the KT theory we make a small digression to discuss the two- 
dimensional Sine-Gordon (SG) model. Its Lagrangean can be written as 

where a is the dimensionless mass parameter, a is a constant of dimension mass~^ 
and P is the SG coupling. It is also integrable and its spectrum and S-matrix was 
also found in 0. The spectrum depends on /3 in a complicated way but it becomes 
simple for the range Stt > > 47r when it is free of any bound states and consists 
of a single 0(2) vector of massive particles whose S-matrix can again be written as 
(0) but now 

ai (6) = cosh smh ■ ' 



2 2 sinh^^^' 

a,(^) = -sinh^sinh^^|i^, (7) 
^ ^ 2 2 cosh^i^ ^ ' 

a^[e) = cosh — cosh ■ ^ 



2 2 cosh^^^ 



where we have introduced the parametrization 

Stt 

The 'isospin 2' phase shift for the SG model is 



sfl{9) = - exp !^2i ^sincj^A;(cu)| (9) 
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with 



• 1 TTUl(l — h') 

~ , ^ smn — T. — - , ^ 

k(uj) = 2v- -. 10 

^ ' 2cosh^sinh|^ ^ ^ 

Note that in the jS"^ 8tt [u ^ 0) hmit the SG S-matrix coincides with the n ^ 2 
hmit of the 0{n) S-matrix, in particular limy^o k{uj) = hm„^2 -^n(i^)- 

The identification of the XY model with the z/ — > limit of the SG model is 
surprising since in this limit the bootstrap S-matrix becomes SU(2) symmetric, 
coinciding with the S-matrix of the SU{2) chiral Gross- Neveu (CGN) model |]TU 



It is not obvious where this enlarged symmetry comes from. The existence of a 
nontrivial XY model is even more surprising in the light of the fact that the beta- 
function of the coupling g"^ in (|l|) vanishes for n = 2 and by making the substitution 

= cos ip , S"^ = simp the Lagrangean (|lD naively becomes free. 

Kosterlitz and Thouless |I| argued that the fact that ip is a. periodic (angular) 
variable plays an important role and therefore the model has nontrivial dynamics. 
They have shown that topologically nontrivial objects, vortices, are present in typical 
spin configurations and their interaction makes the theory nontrivial. 

The standard lattice action of the XY model is 

Sxy = kJ2[^- cos {ipix) - ip{x + A))] . (11) 

X,fJ, 

We denoted by K the inverse of the XY model coupling to avoid confusion with the 
SG coupling j3. Assuming universality, not only the cosine function but any other 
27r-periodic function W{p>) which has a local minimum at = defines a possible 



XY model lattice action. The Villain model action [|TT| is chartacterized by 

W{ip) = -^ln[j2ew{-^iip- 2vrm)^}] . (12) 

Kosterlitz and Thouless showed that typical spin configurations can be repre- 
sented as a mixture of smooth, topologically trivial configurations (spin waves) and 
a gas of vortices (of integer topological charge). The KT vortices are not interact- 
ing with the spin waves, but there is a logarithmic interaction potential between 
the vortices which are therefore identical to a two-dimensional Coulomb gas. This 
spin wave + Coulomb gas (SWCG) picture is only approximate if we start from the 
standard action ([TTD but it is an exact duality transformation for the Villain 



action corresponding to {^)- That the XY model with standard action is in the 
same universality class as the Villain model was demonstrated using Monte Carlo 
renormalization group techniques [|T3|. On the other hand, it has been shown rigor- 



ously [|l^ that the Coulomb gas has a phase transition point at some finite critical 
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coupling Kc- KT interpreted this phase transition as one of vortex condensation and 
by a (heuristic) energy-entropy consideration showed that in the vicinity of vor- 
tices of topological charge ±1 only are important, higher vortices can be neglected. 
It is easy to see that this system (SWCG with unit charge vortices only) is exactly 
equivalent to the SG model. In ref. it was shown that the extremal SG fixed 
point l3* = v^87r,Q;* = is appropriate to discribe the KT phase transition. The 
renormalizabilty of the SG model around this point was explicitly demonstrated up 
to two-loop order in a simultaneous perturbative expansion in a and S = ^ . 

Finally, there is a further transformation that explains the dynamical SU{2) 
symmetry of the XY model. The SG model can be exactly mapped to a fermionic 
model formulated in terms of a two-component Dirac fermion ip. The transformation 
is similar to the well-known one that relates the SG model to the massive Thirring 



model [161. Here the fermionic model is a deformation of the chiral Gross- Neveu 



model with four-fermion interaction: 

t1\2 „ I i2\2 ( „ , £ \/ t3\2 



Cf = ^(^7AV^) - 9o{r,r - 9o{4r - (90 + h){4Y , (13) 



where 



= I ^7m^> (14) 

is the fermionic SU{2) current. The relation between the SG couplings 6 , a and the 
fermion couplings go , /o is 

8 1 

« = -fl'oH , S = (5-0 + /o)H , (15) 

vr TT 

where the dots indicate that the relations (|15D receive higher order corrections in 
perturbation theory. In the fermionic formulation the KT fixed point is the Gaussian 
one and for vanishing deformation parameter, /q = 0, the model is manifestly SU{2) 
symmetric. The corresponding relation in the SG language is a + 85 = at lowest 
order. 

To summarize, the XY model in the vicinity of the Kosterlitz-Thouless transition 
point is believed to be described by the SG model with extremal coupling /3 = v^Svr. 
This is further equivalent to the two-component chiral Gross-Neveu model around 
its Gaussian point. We will use the SG language throughout this paper. 



3 The SG description of the 0(2) model 

In this section we review the SG description of the KT theory closely following the 
approach of Amit et al. 0]. Without loss of generality we can adopt the somewhat 
unusual regularization scheme of the authors, since, as we will see, all important 
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results are universal, i.e. independent of the regularization scheme. Nevertheless, 
it would be interesting to repeat all the calculations below using some of the more 
customary regularizations like the lattice or dimensional regularization. 
Our starting point is the Euklidean Lagrangian ^ 

L = ^d,<pd,<P + ^ ^ [1 - cos(/3o0)] , (16) 

where tuq is an IR regulator mass and a is the UV cutoff (of dimension length). We 
have denoted the dimensionless SG couplings by /?o and to emphasize that they 
are bare (unrenormalized) . UV regularized correlation functions are calculated by 
using 

Go(x) = ^Ko (mov/^^T^) (17) 

where Kq is the modified Bessel function, as the propagator. Our strategy is 
slightly different from 0, who really considered the renormalization of the massive 
SG model ( [TB| ) of mass mo. We treat rriQ as an IR regulator mass and consider IR 
stable physical quantities for which we can take the limit tuq — > already at the UV 
regularized level (before UV renormalization). All renormalization constants are, 
for example, IR stable and independent of ruQ. 

The SG coupling Po is close to its special value \/87r and a simultaneous pertur- 
bative expansion in the mass parameter and the deviation 6o is defined, where 

/52 = 87r(l + 5o) (18) 
and the parameters are renormalized according to 

ao = ZaCK , (19) 
/5o' = Z-'(3\ (20) 

Here the Z-factors are functions of the renormalized couplings a and 5 and the 
combination 

I = In fia, (21) 

where /i is an arbitrary mass parameter (basically the normalization point). Sim- 
ilarly, a renormalization constant Z is necessary to make Q, the spin-spin 2-point 
function finite: 

Gb = Zg . (22) 

For vanishing mass parameter a the Lagrangian (|16D is trivial and Z^ = 1 since 
there is no need to renormalize the SG coupling. The spin-spin correlation function 



6 



(which is an exponential of the basic field (p) gets renormalized even in this point, 
but in this case its renormalization constant is simply 

Z = (/ia)~OT) . (23) 

In addition, there is a symmetry a —a (which corresponds to a ^ shift in the 
basic field). 

Taking into account the above constraints, the perturbative expansion of the 
Z-factors must be of the form 

= l + f,a'l + a^6{f2f + f2l) + ---, (24) 
Z„ = l + g,6l + a'' (^2/' + 92l) + 5' (^3/' + 93l) + -- - , (25) 
Z = e~^{l + a^{hil^ + hil) + ---} . (26) 

Amit et al. found the following results. 

f = - f = f = -A 

1 5 
^i = -2, 92 = ^, ^^"~64' ^'^^^ 

93 = 2 , 5,3 = . 

Furthermore hi = —1/256, but the number hi is not known at present. The above 
two-loop beta-function coefficients were calculated also by other methods. The 
original calculation has recently been reconsidered and the results (^) have been 



confirmed [17 



The spin-spin correlation function satisfies the equation 

vg = jg, (28) 

where V is the renormalization group (RG) operator 

d Q Q d 

^ " '""d^ ^ ^ " """^^ ^ ^° W ^ ^^^^ 

and the /3 and 7-functions are given by 

/3a(ao, 5o) = -9iao5o - 92al - S'sao^o ' (30) 

l3s{aQ, 5q) = fial + (/i + /2)ao^o H , (31) 

7(«o,^o) = -^ + ^'^o-^'^o + ^i«o + --- • (32) 

Now, as is well known, not all /9-function coefficients are universal. For example, 
under a redefinition 

ao = ao + LaoSo + ■ • • , (33) 
~5o = 5o + Kal + --- , (34) 
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the coefficients change as 

91 = 9i, fi = fi, 

h = g2-Kg,-Lf^, h = h-2Kg^-2Lh, (35) 

h = 93- 

(( ^3| - ^ ) is the most general perturbative redefinition respecting the ao ^ —ao 
symmetry together with the requirement that for = (5o is not redefined.) From 
(|35|) we see that in addition to the one-loop coefficients gi and /i there exist also 
two-loop invariants. They are and the combination 

J = 2(72 - /2 . (36) 
Two important physical quantities are the correlation length 

^ = -L = e*{°o,5o) ^37) 
^ Ma ' ^ ^ 

where M is the mass of the physical particle and the dimensionless susceptibility 

X = ^ f rf2xe?(x) = e*(°0'^o). (38) 
J 

From (pSD it follows that the exponents satisfy 

= 1, (39) 
D$ = 2 + 7. (40) 

It is useful to introduce the RG invariant quantity Q{ao, 60) which satisfies DQ = 
0. Introducing the inverse function k{6Q, Q) that satisfies 

ao = k{6o,Qiao,6o)) , (41) 

we can define new (3- and 7-functions: 

b{So,Q) = (3s{k{6o,Q),So), (42) 
r(5o,Q) = ^s{k{So,Q),5o). (43) 

The advantage of using the variables 60 and Q is that the RG invariant Q can be 
treated in many respects as if it were a numerical constant and were a single 
coupling constant. For example, if we write 

^(ao, So) = ^i((5o, Qiao, So)) (44) 
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and 

<l>(ao, 5o) = $i(5o, <5(tto, 5o)) (45) 
then the functions and $i can be determined from 



and 



^'i'i(c>o,Q) = ^7^-7^^ (47) 
o6o b{6o,Q) 



respectively. 

Using (pOD and ([3T|) Q can be determined. 



Q(ao, 5o) = - 25o + 2^72ao'^o + F26I + ■■■ (48) 



and using this in (|42|) we find 

K<^o, Q) = Q + 25o' + AQ5o + B5I + --- , (49) 
while the F function to this order is 

r(5o,Q) = -7 + 7^o + ---. (50) 



Here 



4 4 



J 2 2 2 4 

^ = 1--^, B = --g, + -A, F2 = -gs + -A (51) 



or numerically 



A = 3, B = 2, F2 = 4, (52) 



but we will keep these constants in the following to explicitly demonstrate that all 
our results are universal. 

We know that a one-parameter renormalizable subspace in the Sq-Oq plane is 
equivalent to the SU(2) chiral Gross-Neveu model. (This is most evident in the 
fermionic formulation.) This subspace must correspond to the Q = RG trajectory 
because we know that it goes through the point ao = Sq = 0. Moreover, it must be 
the (5o < branch of the Q = trajectory, since it is the one that is asymptotically 
free in perturbation theory. 
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Figure 1: Phase diagram of the KT theory. The entire Region I is criticaL The model is massive 
in Region III {Q > 0) and Region II {Q < 0). The dotted lines are RG trajectories. The full 
lines and S2 correspond to the critical surface of the XY model and the (bare) CGN model 
respectively. The dashed curve is the XY model hitting the critical surface at c. 

Following IP the phase diagram of our model is represented in Figure 1. The 
CGN model corresponds to the separatrix S2 on this plot, which will be referred 
to as PD for short. Region III corresponds to Q > 0, whereas Regions I and II 
correspond to Q < 0,6o > and Q < 0,6o < 0, respectively. 

In the neighbourhood of S2, i.e. close to the CGN case we have 

b{So,Q) = bo{6o) + Qh{6o) + -- - , (53) 

where 

60 (5o) = 261 + B61 + ---, (54) 

h{5o) = l + A6o + ---. (55) 

It is a crucial observation that near the CGN line the correlation length exponent 
is a smooth analytic function of Q: 

^i((5o, Q) = Ho{6o) + QH,{So) + ■■■ , (56) 
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where 



HoiSo) = -^-^\n\26o\ + ^o + --- , (57) 

We will calculate the value of the nonperturbative constant \l/o in Section 5. 

Now we can integrate p6|) using the perturbative expansion p9|) and for Q > 

we get 

^,{5o,Q) = ^( l + arctg^j -|ln2(Q + 25o') 

Q_ 

8 g + 26^ 




+ o ^ .,2 +^o + ---- (59) 



Here the dots stand for higher order terms in the perturbative expansion (these are, 
in principle, calculable) and also for an unknown (nonperturbative) function of Q 
only. For small Q (|59|) becomes 



^,{6o, Q) = ^ 0(<^o) + Ho{5o) + QH,{5o) + ■■■ , (60) 

where the dots stand for terms higher order in Q. They come from the higher 
perturbative terms of ( |59D and also from the nonperturbative function mentioned 
above. The point is that there are no terms, singular in Q, coming from any of 
these two sources. This is obvious for the perturbative terms, but must also be true 
for the nonperturbative contributions, since otherwise would be singular on the 
5*2 line. Requiring it to be nonsingular on 5*2, we force to diverge on 5*1, which 
is therefore part of the critical surface of the phase diagram PD. An other line, on 
which we know the correlation length must diverge is the 6q axis ao = 0, because 
this corresponds to a free, massless model. 

For Q < 0, it is convenient to parametrize Q in terms of the 5q value at which 
the RG trajectory intersects this axis. In other words, we have to express Q in terms 
of d that solves 

b{d,Q) = 0. (61) 

(Note that in this parametrization \So\ > \d\, because of Og > 0.) The perturbative 
solution of (|6lD is 

Q = -2d^ + {2A - B)d^ + ■■■ . (62) 
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Using this parametrization the perturbative solution for Q < is 



^i(^o,Q) = -^ln( 



I . ^6o + d\[ 2A-3B 1 2A-B d 

; 1 + d + 

6o-dJ I 4 J 8 6o + d 

j\n2\6o + d\+y{d) + --- , (63) 



where 



for c/ < y{d) = + ■ • • , (64) 
for c/ > y{d) = oo . (65) 

(|6^) is obtained by matching (^) to (|56| ) for small Q, while (|65D is formally true since 
this is the only way to achieve that the correlation length diverges on the positive 
half of the 6o axis. (Without the infinite constant ^ actually vanishes there.) 

Now we can discuss the phase diagram of our model. The entire Region I is 
critical. The massive phase is Regions II and III and the critical surface bordering 
them is 5*1 plus the negative part of the 6q axis. They are smoothly connected 
across 5*2, which is the (bare) CGN model. The 0(2) NLS model corresponds to 
the dashed curve of PD. In a MC experiment we are approaching the critical point 
c from the massive phase (Region III). We will denote the 6q coordinate of c by 
do- Because the RG trajectories are running basically parallel to Si, it is physically 
irrelevant at which point the critical surface is reached and therefore the parameter 
do is irrelevant. The continuum model will be the same for all points on 5*1, including 
the origin. But the origin is the point, where (coming along S2) the continuum CGN 
model is defined! So our continuum theory is inevitably identical to the (massive 
part of the) S'?7(2)-invariant CGN model. 

If we start from somewhere in the middle of Region II, we can define a massive 
continuum limit by approaching the negative half of the 60 axis. The intercept d is 
then relevant. The continuum theory is the SG model with 

P^ = 87i{l + d). (66) 

Returning to the 0(2) model, the dashed trajectory can be parametrized as 

So = do + diT H , (67) 

ao = k{do,0) + aiT + ■ ■ ■ , (68) 



where 



K,-K (69) 
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is the reduced coupling and we have assumed that physical quantities are analytic 
in K. {Kc = 1.1197(5) JlSl.) Then also Q is analytic in r: 



g ~ r . (70) 
From (|60D we see that along the 0(2) curve 

lne=-|= + ifo(rfo) + ---, (71) 

where the dots stand for terms analytic in Q (and vanishing for Q = 0). In other 
words, for the 0(2) model [|l| 



^ = Oexp{-^|(i + o(^)), (72) 



where the constants C and b are not universal. This is the famous KT formula 
showing that the phase transition is of infinite order in the reduced temperature. 
It is more important for us that ([7T|) can be rewritten as 

2(ln^ + u 

where 



g = , + (73) 



u = -Ho{do), (74) 
which is given (if do is sufficiently small) perturbatively by 

1 B 

Trr + -rH'^do)-^o. (75) 

Note that the leading l/(log.^)^ term in (|73D is universal and only the subleading 
terms (depending on the value of the parameter u) are model-dependent. 

The susceptibility exponent $i can be studied similarly. It can be written as 

$1 = I + $1 = I VI/, + 1 ln(252 + Q) + . . . , (76) 



where $i satisfies 

d^i _ T{6o,Q) + l 
dSo b{6o,Q) 

Now the crucial observation is again that, for small Q, 



(77) 



7 

$1 = - ^1 + ci + csQ + ■ ■ ■ , (7J 
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because the (calculable) perturbative terms are analytic in Q, while the (not calcu- 
lable) purely Q-dependent terms in $i must also be analytic in Q otherwise these 
latter singularities would also turn up for the CGN line 5*2, where they must not. 
From (ffHl) we have 



X 



e^e^^(l + C2Q + ---), (79) 



i.e. there are no (multiplicative) log corrections in the scaling relation for the sus- 
ceptibility. The possibility of such multiplicative logarithmic corrections is discussed 
Pi. 



m 



4 Determination of the lattice artifacts 

Recall that the RG invariant Q has a completely different meaning for Region III 
(which contains the massive phase of the XY model) and for Region II (where the 
usual massive SG model with < Stt can be defined). Indeed, in the positive 
(5o part of Region III, close to Si, the (positive) parameter Q merely measures the 
distance from the XY critical surface on which it vanishes, whereas in Region II Q is 
a relevant (negative) parameter related to the SG coupling (3 by (|6T|) and (|66|). Our 
main assumption is that in spite of this difference physical quantities are smoothly 
depending on Q in the vicinity of the separatrix 5*2 connecting the two regions. More 
precisely, we will assume that close to S2 any physical quantity U has the form 

U{Q) = Uo + UiQ + 0{Q^). (80) 

Here Uq = U{0) is its value for the CGN model (and thus also in the continuum 
limit of the XY model). The first correction coefficient Ui can be calculated from 
the SG model as follows. Using the identification ( |6TD and its perturbative solution 
(|6^) together with (|66D and (j|) we have 



g = -2z/2 + 0(z/4). (81) 
This means that if in the SG model, close to the CGN point z/ = 0, we have 

U{u) = uo + uiu^ + 0{u''), (82) 

then 

Uo = uo and Ui = -ui/2. (83) 

Translated to the language of lattice artifacts by ([73|) we thus have 
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This means that lattice artifacts typically go away very slowly, only as l/(log.^)^. 
On the other hand the leading artifacts are universal and calculable. 

We apply this method first to the scattering phase shifts. Recall the SG model 
S-matrix (0) with (0). The three distinct S-matrix eigenvalues are 

4?(9) = 2^,(9) + = i^M(!l±||4=)(9,, (85) 



and 



.(2) 
'SG 



(e) = (T2ie) + a^ie) = - exp |2z ^ sm{euj) A;(cj)} 



^7) 



where k{u) is given by (|TUp. 
We now write 

s« (^) = exp {2z5jV) + 2iu^S?{e) + 0{u^)} (88) 

for i = 0,1 and 2. Here ^q*"* are the CGN phase shifts which, as remarked before, co- 
incide with the n — i> 2 limit of the 0{n) phase shifts. The first correction coefficients 
can be obtained by a simple calculation. The result is 

<>W = 0. = = (89) 

This can be used to obtain the leading lattice artifacts in the XY model by the 
relation (18^ 



5 Current-current 2-point function and free en- 
ergy 

Consider the 2-point function of the Noether current 

Jf, = i^e^^d^(f) . (90) 
Ztt 

Its Fourier transform I{p) is defined by 

{Mx)My)) = 1^2 e--(-^)^^^^^i^^ i{p) . (91) 
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It is easy to calculate I{p) to second order: 

I{p) = ^l^l + 5o + al[nip,a) + --- }, (92) 
where 

n{p,a) = Mlnpa + no) (93) 

and 

Qo = C-^-\n2, (94) 

C being Euler's constant. 

Standard RG considerations give 

I{p, 5o, Q, a) = I{po, A, g, a) = i{p/M, Q) , (95) 

where the running coupling is the solution of 

^i(A,Q) = ln^-lnpoa. (96) 

Let us consider the Q < 0, 5o < case (Region 11) first. For p ^ oo also \t'i ^ oo 
and therefore A — > rf, where d is defined by 

h{d, Q) = or equivalently ao = k{d, Q) = . (97) 

This gives 

/(oo,Q) = -(1 + d), (98) 
7r 



where d parametrizes Q according to (|^)- (EH) is consistent with the identification 

Next we study the case of small Q, because this is relevant if we are interested 
in the approach to the continuum limit along the 0(2) curve. We assume that / is 
analytic in Q: 

i{z,Q) = io{z) + Qh{z) + --- . (99) 

It is a standard exercise to obtain the asymptotic expansion of the coefficients Jq 
and Ji in perturbation theory. One gets 

Jo(p/M) = ^|l-A + 2/€A' + ---|, (100) 

~ , , 2 f 1 2B-3A , , , 
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where 

and the couphng A is the solution of 



(102) 



(103) 



The asymptotic expansions ( |100[ ) and ( |101D are vahd for A -C 1, i.e. for p — » oo but 
also Q must satisfy 



g<6A2 



(104) 



so that the expansion ( p9D makes sense. 

It is by now standard how the nonperturbative constant k can be calculated. For 
this it is necessary to consider the free energy in an external field and we now turn 
to this calculation. We follow here |19] and start from the modified Lagrangian 



Ch = ^{d^(t)f - cos(/5o0) + ^<9i0 , 
2 PgO^ 27r 



(105) 



which corresponds to adding a term ihJ2 to the Lagrangian density. The modified 
ground state energy must be of the form 



(106) 



where J-'{h) is dimensionless. 
In perturbation theory we get 



where 



and 



J^{h) = l + So + al n(h,a) + 



Cl{h, a) = /i(ln ha + CIq) 



fto = C -- + ln2. 



(107) 



(108) 



(109) 



For Q = therefore J^{h) has the following asymptotic expansion 

1 5lnln^ K + ln4 



J^ih) = 1 



2 In — 



8 1"^^ 



+ 



2 In 



2 h_ 

M 



+ 



In In — 

lii Hi 



In' 



3 b_ 

M 



;iio) 
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This is to be compared to the exact result ||20| 
T{h) = 1 



1 Inln — 
111111 

2 In — 4 In^ — 

z in ^ 111 j^,^ 



+ 



§ln2-i-iln7r , 



In In — 
111 111 



21n^ — 



111) 



It is a nontrivial check on the overall consistency of our considerations that using 
(|5^) for the numerical value of our universal constant B the third term in (|110|) 
matches the corresponding one in ( |111|) . Comparing ( |110| ) to ( [111|) we also get 



K = In V 27r . 

4 



(112) 



Using the exact value ( p.l2| ) we obtain the asymptotic expansion of the first correction 
in Q: 



Note that the leading correction term (|113|) does not contain any free parameter. 
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